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GROUND LEVEL ACOUSTICAL FOCI IN A 
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By 

Willi H. Heybey 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 

/ 4 /°° 

An infinitesimally slender bundle of sound rays emitted by a point 
source can, under certain meteorological conditions, be returned to the hori- 
zontal plane on which the source sits, and at the same time, may be narrowed 
down to enclose an area actually equal to zero. In such cases the common point 
of arrival is called an acoustical focus. The intensity of returned energy there 
attains a very high degree capable of damaging effects. 

For determining such points the atmosphere is usually divided up into 
layers in each of which the change of the sound speed with height is considered 
constant. Three layers often suffice to approximate the many layers arising 
from observation. 

The report studies focus formation engendered by the third layer in 
collaboration with the two lower layers. A focus may appear in variegated cir- 
cumstances. For each case a computational scheme is given which may be 
followed by slide rule or assigned to machines if higher accuracy or a system- 
atic survey is desired. 
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TECHNICAL MEMORANDUM X-53344 


GROUND LEVEL ACOUSTICAL FOCI IN A 
THREE -LAYERED ATMOSPHERE 

SUMMARY 


The report classifies the states that may occur in a three-layered at- 
mosphere. 

For each such class the possible existence of foci is discussed in 
developing the iterative method for solving the equation valid for foci that evolve 
from the third layer. Those originating with the second layer are not considered; 
for these, existence conditions and modes of computation are simple and known. 

The location of the focus is given, and the seriousness of the intensity 
near it is described by a characteristic quantity. 

Unexpected features of focus formation are pointed out as, e. g. , exis- 
tence of as many as two third-layer foci in conditions where a similarly struc- 
tured two-layer atmosphere would support none (even otherwise, it could produce 
at most one) . 

Three classes permit a non-iterative, direct solution of the focal equa- 
tion. They throw in vivid relief some of the novel aspects including the way in 
which the focal distance and intensity respond to variations of the third-layer 
velocity gradient. 


I. INTRODUCTION 


The complex propagation patterns of sound energy in the natural atmos- 
phere become amenable to relatively straightforward mathematical description 
when three major simplifications are made. First, of all the effects a sound 
wave may experience, refraction and reflection only are admitted ( ray acoustics) . 
Secondly, the propagation velocity of sound is considered a function of height 
only ( stratified atmosphere) . An acoustical ray sent out by the point source in 
a given azimuthal direction then will forever remain in the vertical half -plane 



fixed by that direction so that the ray pattern in it can be determined without re- 
gard to other such planes. Thirdly, the wind must not have a vertical component 
(which, incidentally, is hard to measure) . 

Rothwell (Ref. 1) has dealt with the problem on this basis. The mathe- 
matical formulation is still rather involved and the final results are not essentially 
different from those obtained when two more restrictions of lesser import are 
made: on the one hand, the wind velocity is assumed as small when compared to 
the thermodynamic sound speed; on the other, the sound ray inclination should 
not go above, say, 6 = 18. 5° (cos 0 = 0. 95). These rays account for the energy 
return to be expected in the nearer vicinity of the source. 

Meteorological observations of temperature, wind velocity, and humidity 
at specified heights can then be compounded to determine the local velocity of 
sound prevailing in a selected half -plane. It is customary to connect the dis- 
crete data points by straight line segments and thus to break up the stratified 
atmosphere into plane-bound layers in which the velocity, V, of sound is a linear 
function of the height (y) . In the kth layer it may be written as 


V=V k-l + 


M k (y 


r k-i' 


(i> 


where the index (k-1) refers to the bottom of that planar layer. The coefficient 
ju k is evidently given as 


M, = 

k y. 


V. -V. , 
k k-1 


k y k-i 


(2) 


and thus is determined by the top and bottom speeds and the layer's vertical 
extension. Its numerical value can be positive, zero, or negative, but is never 
-2 1 

large in practice (10 in general) . The sound source is placed at the height 

S6C 

y 0 = 0, so that in the ground layer 


Mi = 


Vj- Z- Yn 

yi 


The velocity polygon (customarily height is plotted versus velocity) is 
different for each vertical half-plane through the source (except in a still 


2 


1 


* atmosphere) . For energy considerations it is necessary to introduce an adjacent 
half -plane in which the state of the atmosphere is substantially unchanged. The 
sound field within the infinitesimal vertical wedge formed in this way is con- 
veniently divided into ray tubes originating at the source and widening or narrow- 
ing on their courses. Through all the rectangular cross sections of a tube, the 
same amount of energy is transmitted per time unit because dissipation is not 
taken into account in ray acoustics. If the cross-sectional area shrinks to zero, 
as it does when the top and bottom surfaces of the tube gradually approach each 
other and finally intersect at zero angle, the intensity (energy/ time area) be- 
comes increasingly larger and at the end infinitely large. In die planar 
representation one then obtains what is called an acoustical focal point. 

Infinite intensity will hardly be observed in practice because, for one 
reason, one cannot measure at a mathematical point; for another, the irregu- 
larities of the natural atmosphere tend to dislocate the focus, if not to obliterate 
exact zero-angle intersections altogether. Also, interference processes might 
weaken the energy return. However, if ray acoustics indicate a focal point the 
prima facie assumption appears justified that high sound pressure levels may 
be expected to occur there. 

Evidently a continuous sequence of zero-angle intersections can exist 
in the open atmosphere forming a spatial point array in the nature of a ray 
envelope. These focal curves are of special interest in the air space over 
slightly hilly terrain and will be subject matter of a later memorandum. The 
present investigation is concerned with their zero-height points only, hence- 
forth called "the" focal points. 

The velocity profile as it is obtained from meteorological observation 
is a many-sided polygon corresponding to a multi-layered atmosphere. The 
raw material is corrected when necessary (cleaned of doubtful points, etc. ) 
and often smoothed (for a more even data distribution) . It may include local 
effects not characteristic for the atmosphere as a whole, and it is subject to 
an error margin as any measurement is. The "observed" velocity profile must 
therefore be envisaged as a more or less accurate approximation to the true 
vertical velocity distribution. 

Since the necessary computer evaluation is numerical, conclusions of 
a general nature can seldom be made. In freak cases, the result may even be 
misleading; an instance is known where, with the usual step increase of A0 g = 

0. 20, a focus, known otherwise to be present, was not detected before the step 
was drastically decreased by the factor 200. 
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Experience has also shown that ray pattern and focus formation can be 
surprisingly sensitive to small variations in the meteorological parameters. 

The first indication of this was given in Reference 2. Reference 3 studies the 
effect of small data modifications in the two-layer set-up where closed expres- 
sions are available eliminating the need for numerical specimen calculations. 

The present report undertakes to open the way for similar investigations 
in the three-layered atmosphere, as far as focus formation is concerned; sys- 
tematic surveys will be made subsequently. Moreover, a three-sided velocity 
polygon can frequently be substituted for the many-sided polygon of observation, 
without stretching the evidence beyond the breaking point. It then may help to 
show that an observed focus could indeed have been present, although the com- 
puter deck did not, and perhaps could not, reveal it. Finally, an insight in 
characteristic features of focus formation will be provided that goes beyond 
what is known from the study of the two-layered atmosphere. 


H. BACKGROUND AND GENERAL INFORMATION 


When ascending rays penetrate into a layer, their upward swing is not 
being checked unless the velocity gradient, m> there is positive. If the layer is 
sufficiently extended in height, the ascent of low inclination rays may be stopped 
altogether and converted into a downward trend; if so, these rays are sure to 
hit the ground again. Among them, there might be a focal ray, that is, a ray 
thatat zero level intersects with the adjacent ray under zero angle. Ray inclinations 
will be identified here by the angle, 0, they make with the local horizontal line. 
Along the descending portion of the path 0 is negative. A ray may be singled 
out of the mass of all rays by its angle of departure at the source, 0 = 6 0 . That 
of the focal ray will be called 0 0 ' r . 

In order to have a focus out of the third layer, the fundamental require- 
ment must be made that 

M 3 ^0, (3) 


which, however, is not the only one. Sufficient conditions will be set up as if the 
third planar layer extended to unlimited heights, allowing all the rays that enter 
it to return to ground level. No conceivable focus then can escape attention. 

After ascertaining that a focus, i.c. , a real value of 0 O * , can exist with a given 
M 3 , one can determine without difficulty whether or not the focus will in fact appear 
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if the third layer is terminated at y = y 3 . The top velocity, V 3 , follows from 
expression (2) . Unless 

cos0 o *£^, (4) 

v 3 

the presumptive focal ray will not be turned around within the third layer and 
will enter the fourth, so that the focus cannot materialize. 

The ground layer can never induce focus formation, but the second 
(middle) layer can, granted sufficient thickness and a positive velocity gradient. 
Second layer foci can be dismissed here, since Reference 4 deals with them in 
considerable detail. The same report also contains the fundamentals of third- 
layer relations, especially those numbered ( 55) , (57), (58), (59) , (p. 40-41 ) . 
They are given in terms of the quantity 

£, = cos 2 0 O ‘' . 

However, it has been found that the mathematical discussion is facili- 
tated if, instead of £,, the cotg - function is used: 


q = cotg 2 0 O * 


(5) 


Equation (55) of Reference 4 contains a misprint (the factor K 2 2 in the third 
square root should be K t 2 ) ; after correction is made, it assumes the form 


1 - 


Do 


\l 1 + 6^q 


D? ~ Da 
\/ 1 + <5 2 q 


= 0 


where 


( 6 ) 


D. 

l 



6. 

i 



K. 

i 


V. 

l 

V 0 




(6a) 
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If a real and positive solution q exists to equation (6) , the third layer induces a 

focal point, the focal ray departing at the source with the elevation 0 O * = arc 

cotg \T q~ . This ray impinges on the horizontal plane at a distance, x , from 

F 

the source given by expression (57) , now written as 


x 


F 


_gyi J_ 

K t - 1 \Tq 




1 - D2 \! 1 + 6jq + (D2 - D3) \TT~+~ 6 ^q 




(7a) 


It is sometimes convenient to use a second representation of x^ obtainable by 
combining relations (6) and (7a) to give 


x 


F 


2 Yi 

Kj ~ 1 


\Tq~ 


~D? d 1 + (D? ~ Da) &> 

\/ 1 + djq s] 1 + ^q 


(7b) 


A measure of how strongly the intensity becomes infinite at the focus is 
provided by the quantity, o> , which is proportional to the intensity average in 
a differential vicinity of the focus; relation (58) of Reference 4 is transformed 
into 


U) = 


4 n/ q 
d 2 x 

X F S 


d ft 


0 


On — On' 


( 8 ) 


where, by expression ( 59) , 


d 2 x„ 


ddn 


On — On 


= . -gii— c 1 + q ) 2 

n/ q 


Ki - i 


+ I P 2 ~ D s) § 2 
(1+ 62 q) 3 / 2 


~ D ? . 

(1 + 6 iq ) 3 / 2 


(9) 


Formulas (7b) , ( 8) , (9) indicate that a; is mainly determined by the 
\l Q 5 ;< 

ratio ( i +"q)'^ * Before e o waxes larger than permissible, the severity of 
a focal situation will in general increase as Oq increases; larger elevation 


6 


focal rays threaten higher sound pressure levels. There are exceptions to this 
rule, as will be seen in Section VII. The focal point generally will move toward 
the sound source if 6 0 ' fi is increased. 


m. ON COMPUTING SOURCE LEVEL FOCI EVOLVING FROM 
THIRD-LAYER CONDITIONS 


In Figure 1 some forms of the three-sided velocity polygon are sketched 
(many more of them must be distinguished) . The slope, tan (ZL - arctan ^ ) , 

2 k 



FIGURE 1. SOME FORMS OF THREE-LAYER VELOCITY PROFILES 
of the kth side is positive, infinite, or negative depending on whether n, is 

K 

positive, zero, or negative. By condition ( 3) that of the uppermost side must 
not be negative. Otherwise, the polygons can have any shape as the choice of 
Vo, Vi, V 2 , y lt y 2 is completely free, at least mathematically. In physical 
reality, the three velocities will not differ very much from each other so that 
by definition (2) , the n as a rule will be small quantities. Likewise, the 

t\ \ 2 

<5^ = 1 - I ~ — I will be small . 
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It is seen that the quantities 6 t and <5 2 can be chosen independently from 
each other, and also from both D 2 and D 3 : 



Note that M3 cannot be replaced here by — ^~ 2_ , since the third layer has been 

Y3 “ Y2 

taken as infinitely extended; the physical impossibility of this cannot prevent us 
from ascribing a value to M3. This quantity must be taken as a parameter in- 
dependent of the five cited above. 

Since the focal equation (6) contains four independent parameters, the 
general discussion of possible roots q = 0 offers great difficulties. It was de- 
cided not to remove the irrationalities by two times squaring: First, the result- 
ing rational equation is of fourth degree so that the roots depend in a complex 
manner on the coefficients; second, the coefficients in themselves are unwieldy 
power combinations of 6 2 , D 2 , (D 2 - D 3 ) , all but preventing an elementary 
classification of the solutions; third, the roots would have to be checked back 
with the original equation because some of them will not solve it. "' 

If one stays with the original equation, the convergence criterion of 
Newton's approximation often requires here that the root sought should already 
be almost exactly known. By contrast, the convergence requirements when 
solving iteratively proved to be less stringent. 

In order that a process 


should converge, it is necessary and sufficient that the root and all the iterated 

Graeffe's numerical method has been successfully programmed for the fourth 
degree equation by Mr. Q. Peasley. It gives the four roots with great accuracy, 
unless two of them have equal modulus, when the deck is not applicable. 
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values leading to it stay in an x-interval where 


V 


If one attempts to satisfy this condition, easier criteria are obtained when the 
convergence is monotonic rather than oscillatory. To insure this, the deri- 
vative must not be negative within the interval. Depending on the signs of the 
coefficients one will set the equation (6) into several iterative forms such that, 
if convergence takes place, it is monotonic. In this context it should be noted 
that the quantities 


u = N/l+6 1 q,v=N/l+6 2 q, r = ^ , s = ^ (1°) 

are all not negative for a physical reason. The ascending focal ray leaving the 
source at the positive ( or zero) angle 0 O ‘ will enter the second and third layers 
at positive angles, 9{'~ and 0 2 ‘ • Physically, the square roots u and v are the 
ratios 

sin sin On " 

sin 6{'' ’ sin 0 2 " 


and therefore not negative. This happy circumstance lessens the variety of 
forms equation (6) must be given and is otherwise of help in solving it. Since 
we wish to stay away from squaring the equation, one of the four quantities ( 10) 
will have to be taken as the unknown variable. From the outset then, the inter- 
val in which to investigate the derivative is restricted to values of the argument 
that are not negative. 

Four basic forms emerge on arranging equation (6) for iterative pur- 
poses: 


r 


_L_ 

d 2 


1 + 


(D 9 . ~ D a) r 



+ 


<5i 


= f(r) 


( 11 ) 


9 
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Many alternative representations can be tried, of which one was found 
to be more useful in certain circumstances than any of the basic forms: 



The roots in these expressions are all not negative, so that the sign of the deri- 
vatives solely depends on those of the parameter combinations in front of the 
roots. Even if these are not negative, the derivative, in addition, must remain 
smaller than unity to insure monotonic convergence. The r- and s- approaches 
are both necessary for this reason, although they have non-negative derivatives 

for the same values of t 2- , D 2 , (D 2 - D 3 ) . 

6 i 

For determining the range of applicability of the several representations 
the quantity 


M2 

has been proved a fitting guide. It can vary from (Mi > 0, ix 2 = 0 ) to +°° 

(/Hj < 0, /u 2 = 0) . This vast interval can be split into subintervals marked off by 
certain values of D 2 for which equation ( 6 ) assumes special forms to be treated 
separately in later sections. Among these values are D 2 = ± » (zero velocity 
gradient in the second layer) , D 2 = 0 (equal gradients in the two lower layers) , 

D 2 = 1 (zero gradient in the first layer; = 0) . To these will be added the value 
of D 2 corresponding to 63 = 0 since this clearly is an exceptional case, too. 

K 2 = 1 here, so that, from the definition (2) 
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This value of D 2 is larger than unity. 


Four subdomains emerge in that way: 


A. 

- 00 < D 2 < 0 

corresponding to 

+00 > to- > 
M2 

1 

B. 

0 < d 2 < 1 

corresponding to 

O 

A 

d £ 

A 

■H 


C. 

1 < d 2 < ^ 
2 yi 

corresponding to 

1 

tH 

A 

£l £ 

A 

O 

Z 2 _ 

yi 

D. 

< Do < +°° 

yi 

corresponding to 

t _Iz > J£l 
yi m 2 

> -CO 


In each of these subdomains the difference 

(D2-D 3 )= (1 -f 2 -) (16) 

M 2 M3 

can have both signs depending on whether M3 > M 2 or M3 < M2- (The case M3 = M2 
is trivial since the atmosphere becomes two-layered.) 

It has been found that in a given subdomain the same iterative repre- 
sentations cran be used to arrive at solutions for q and thus for . A single 
one always suffices when the (non-negative) gradient in the top layer is larger 
than that in the layer below it. One may recall that, in a two-layered atmosphere, 
ground foci can exist on this condition only. It will be shown in Section V that 
with three layers, the condition can be violated without completely eliminating 
focus formation. On the contrary, as an additional novel feature, even two focal 
points may appear in some subdomains of D 2 , so that two different iterative 
processes will have to be used in these. 


IV. VELOCITY POLYGONS WITH M3 > M 2 


The gradient M3 must obey the condition ( 3) ; M2 and m 1 should not be equal 
or zero, nor V 2 equal to Vo, since otherwise equation (6) would assume special 
forms not covered by the representations (11) to ( 15) . 


As an example, the class C will be treated in some detail. Since 

y 1 

< - — — in that class, the ratio 

Mi y 2 - Yi 


&> K, + 1 



ill y? - yi 

Mi yi _ 


(17) 


& 

is negative in it, < 0. Bearing this in mind and applying the approach ( 15) 

°i 

we derive from the convergence criterion g' (v ) < 1, that any iterated value, 
v , must satisfy the inequality v 


(v + D 2 - D 3 ) > - 


61 


Lz_Q_ = M 

Q 1 


) Q 2 (i -Q) 


(18) 


where Q 3 is defined as 


Q 3 = 



(19) 


If all the v do so, the convergence is ensured and will be monotonic since, 
^ dg 

with the approach (15) , is always positive. 

v 


If the positive quantity Q is at least equal to unity, the criterion ( 18) 
is satisfied with any v . However, when 


0 < Q < 1 
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we must specifically make sure that all the iterated values (including the final 
one considered as the solution) reside in the convergence interval (18) . Suppose 
the vth iterated value does so. The (v + l)^ 1 one, 



then has a lower bound 


% + i > ''o-fb o-Q) 

which, to secure convergence, ought to be at least equal to the bound imposed 
by the criterion (18): 

J~( 1 - Ij-Ml -Q) ^-(D 2 - D 3 ) + Q J (1 - ^-) (1 -Q) • 

By expression (16) the difference (D 2 - D3) is positive in the subdomain C as 
long as we keep to the assumption that ju 3 > ju 2 . With Q being smaller than 
unity, the above condition is then satisfied. It would not necessarily be so in 
a case where g 3 < /x 2 which therefore calls for a different iterative approach. 

It remains to determine an appropriate initial value, Vj, to start the 
process when Q < 1. For this purpose we may consider the function 


G(v) - g(v) -v 


which will be zero for the desired solution or solutions of equation ( 15) . The 
curve G(v) has a point of inflection at v = v = - (D 2 - D 3 ) < 0 and a maximum at 


v = v ;,: = - (d 2 - d 3 ) + q ^(i- ^-) ( 1 - Q) 

which, however, does not exist with Q 2 1. Typical forms are shown in Figure 

2 . 
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FIGURE 2. PLOTS OF THE FUNCTION G VERSUS ITS ARGUMENT 
Since in the subdomain C 

d 2 2 ( i - f -) 2 >1 

it follows from the definition ( 1 9) that 



and that, therefore, as Q < 1, 

(1 - |2-) (1 - Q) < 1. 

As a consequence, v* < 1. The convergence condition (18) can be written as 

v > v' 1 ' 

v 

from which relation it appears that a suitable first value could be taken as 

vi = 1 . 


With Q §? 1 the same value will do (as, indeed, any other would) . 
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Figure 2 indicates, in addition, that in all circumstances there will be 
one, and only one, point at which G is zero with a positive value of v. (It is 
easy to back this up by direct analysis of the function G) . Thus, in a three- 
layered atmosphere where p 3 > p 2 and D 2 is located in the subdomain C there is 
always a single focal ray setting out at the source under an angle 0 O * given by 

/ 1 - v 2 

c°tg 9<r 

where v is found by the iterative process (15) . Its point of impact is deter- 
mined by expressions (7a) or (7b). 


Similar derivations can be made with D 2 chosen in the three other sub- 
domains. To summarize: In subdomain 


A. 


B. 


use s-approach and s t = 



use'approach (14) and Vj = 1 + 


D3 - 1 


C. use approach (15) and v t = 1 


D. use r-approach and m = — 2 -~ — A — 

! _ h. D 2 

In all these atmospheric conditions a single zero level focus will exist rooted in 
the (unlimited) third layer. The eight types of velocity profiles involved are 
shown in Figure 3; care has been taken to depict p 2 # p l5 ^ ± 0, n i + 0, V 2 + V 0 . 
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FIGURE 3 . TYPES OF VELOCITY PROFILES WITH M3 > M2 
V. VELOCITY PROFILES WITH 0 < M3 < M2 


There is no chance here that M2 could be zero or M3 negative, but the 
other exemptions listed at the beginning of the preceding section apply also. 

The subdomain C where 


Mi 

M 2 


< 0 , 


^ 2 _ 


< 0 


may again be chosen as a model. Statements more definite than before can now 
be made. Since the second layer gradient M2 is positive, Mi will be negative 
with the further consequences that V t < Vq, Kj < 1, > 0, 62 < 0, K 2 > 1 , V 2 > Vq. 

The polygon C in the later Figure 5 reflects these results. 


Expression ( 16 ) shows that (D 2 - D 3 ) is now negative, while D 2 is 
positive in subdomain C. All of the processes ( 11 ) to ( 15 ) wou ld converge 
monotonic ally. If we select the u-approach ( 13 ), u = \l 1 + <5 jq will be larger 


than unity, and u 2 must be smaller than 



since otherwise the function 


17 


F(u) = f ( u) -u = (Da -u) 


(Da ~ Ds) u 

7(1 _ ' ^ 2 - »* 


( 21 ) 


-) + 


would not be real. If a focus is to exist, a physically relevant solution of the 
equation F = 0 must be sought in the interval 


1 < u < 




<52 


( 22 ) 


Aside from the quantity ( 19) , Q 3 , it is useful here to consider a further 
positive combination, R, defined by 


R 3 = 


(Da -D s) : 

i-A- 


(23) 


The convergence criterion, f'(u )< 1, then requires that 


u < J( 1 - 




) (1 -R) 


(24) 


If u^ obeys that condition, the next iterated value 


( Da - D 3 ) u 


u = D 2 
v+i 


J(i - f") + I 2 " u 2 

Oj Oj V 


is subject to an upper bound 


Vl < D 2 + R (1-R) 

which, to secure convergence, must not be larger than the right side in the cri- 
terion (24) . On substituting Q for this gives the condition 
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Q + R5 1 


(25) 


which must be satisfied by the given values of ■~ 2 - , D 2 , 
here does not take place in all circumstances asht did, 


IV. 


(D 2 - D 3 ) . Convergence 
regarding v, in Section 


Condition (25) has both physical and geometric connotations. 


As to the first, we note that the physical inequality (22) must be ex- 
pected to hold as well at least for iterated values u close to the one that is taken 
as the root u. Now, Q and R, being positive, are true fractions a ccording to 

it is required. Moreover, condition (25) gives 


condition (25), so that, by the convergence criterion (24), u^ < 


as 


(1 -R) (1 - |-) Ir q (1 - &-) 


Since relation (20) is valid in the subdomain C, the right side here is larger than 
unity. Were it otherwise, the left member of the continued inequality (22) 
could not be satisfied with any u^ because of condition (24) . 

Geometrically, analysis of the function F(u) shows that the curve F can 
have forms as depicted in Figure 4. There is a minimum at 


u = <s/( 1 — 


<52 


) (1 -R) 


The value F(u ' ) is positive, zero, or negative depending on whether Q + R = 1. 
It is seen that, if the criterion (25) is satisfied, there exist in general two rele- 
vant roots u of the equation F(u) = 0, except with the sign of equality when there 
is only one. If Q + R > 1, positive roots cannot exist, and focal points will not 
originate in the third layer. It is remarkable that the existence condition for foci 
also insures the convergence of the iterative u-process. However, the latter 
gives only one of the roots. Since the convergence criterion (24) may be written 
as 


u < u ' 

v 
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FIGURE 4. GEOMETRY OF THE FUNCTION F(u) 

this is the one at the left of the minimum abscissa. A suitable starting value 
for the iteration will be the abscissa, u = Uj, at which the tangent to the curve at 
u = 1 intersects with the u-axis: 


D-, - 1 

U 1 = 1 + T (26) 

1 + (D 2 - D 3 ) (1 - ^-) 

Looking for the second root we may try the s-approach ( 12) where the 

convergence criterion, f'(s ) < 1, takes the form 

v 


9 1 Q + R 

s > — — 

v R 1 _ h. 

" <5i 


(27) 


It can be shown that all the iterated values remain in this interval if, again, 
condition (25) is satisfied. This remains true, even when the criterion (27) is 
made stronger by demanding that 
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1 + 6 t q = u 2 > (1- ~ L ) ( 1 - R) = (u*) 2 . 

V V 02 

Thus, the s-approach yields the desired value of q corresponding to the root of 
F(u) = 0 at the right of the abscissa u' ,: . A suitable first value, 8 = 8!, is 
given in the list below. 

The u- and s-approaches converge monotonically, under the conditions 

<5 

stated, not only in C, but also in the subdomain B where, however, ~~ > 1, 

6 i 

and therefore R < 0. 

The study of the velocity profiles with 0 2 /n 3 < n 2 uncovers a less uni- 
form picture of focal possibilities than was found with > /u 2 . A survey follows. 

Subdomain A. (D 2 g 0) : No ground level foci originate in the third layer. 

Subdomains B. and C. (0 < D 2 < ^ , excepting Do * 1) . 

yi 

There is no focus in B, if Q + R < 1 
There is no focus in C, if Q + R > 1 
There is one single focus if Q + R = 1 
There are two foci in B, if Q + R > 1 
There are two foci in C, if Q + R < 1 


i 
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These two can be determined by employing the u- and s-iterations and the initial 
values 


Uj = 1 + 




1 + (D 2 - D 3 ) (1 - % 


Si = 


I>2 ~ D3 


J77J 


- 1 


If Q + R = 1, a direct solution without iteration can be given: 
cotg 2 e 0 ‘ = 


6 , 

2/2 

D2 ( l - 

- 


Subdomain D. (D 2 > ~ There is always one, and only one, focus 

y 1 

originating in the third layer. Use u-approach and 

% 


u i - d 2 - \i ( d 2 - D 3 ) . 


Figure 5 sketches the four types of profiles involved. The type A where 
the polygon sides consistently turn to the left produces no ground layer focus, 
either from the third or from the second layer, although the velocity gradients 
are positive so that all the rays emitted return to the source horizontal. Note 

y? 

also that if D 2 > ( subdomain D) there will always exist one single zero level 

yi 

focus from the unlimited third layer, regardless of all other meteorological data, 
especially of whether M 3 < m 2 or ^ > M 2 (Fig. 3) , -as long as M 3 > of course, 
is not negative. 
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FIGURE 5. THREE-LAYER VELOCITY PROFILES WHEN 0 % < p 2 


VI. THREE SPECIAL TYPES OF VELOCITY PROFILES 

The focal equation (6) for determining 0 O '' assumes different analytic 
forms with those values of D 2 that separate the four subdomains A, B, C, D. 

With D 2 = ±«° and D 2 = 1 a limiting process is called for; with D 2 = 0 and 

D 2 = -^ equation (6) affords simple solutions without the need for iteration. 

The fourth case (D 2 = ^ , V 2 = Vo) is especially instructive, since it comprises 

Yl 

three fairly general classes of three-sided polygons that can be treated system- 
atically without numerical effort ( Section VII) . 

The iteration procedures recommended in the following can be sub- 
stantiated by analyses after the pattern set forth in discussing the subdomain C. 
Special turns have to be added in places. However, as in the preceding sections, 
leaving small conceptual laps open seemed preferable to developing the main 
guiding ideas by often lengthy (and rather elementary) mathematical reasoning. 
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A. D 2 = ± 00 


The gradient in the second layer here is j u 2 = 0, so that V 2 = Vj and 
62 = The profiles associated with these data fall into the two categories in- 
dicated in Figure 6. The case tq - 0 adduces a two-layer situation and is ex- 



FIGURE 6. PROFILES WITH ju 2 = 0 

eluded for that reason. Since both D 2 and 6 2 contain aggregates in K 2 and Kj, 
one performs the necessary limiting process by letting K 2 approach K^. Equa- 
tion (6) then assumes the form 


u 3 = 0 + (D 3 - /3)u 2 

where 

R = K i y ?. zjll 

P K t + 1 y'i 

For iterative treatment one best introduces the unknown 
w — u 3 = (1 + 6jq) 3 / 2 , 


( 28 ) 


(29) 


( 30 ) 


so that 
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(31) 


w .< = /3 + (D 3 -/3) w 2 / 3 = h(w ) . 

Since -y— is positive with (D 3 - j3), the approach (31) will be used only when 

GW 

V 

(D 3 - /3) >0. An initial value has been determined as 
Wi = (D 3 -/3) 3 + 3/3 . 

The very special case (D 3 -/3) = 0 gives the direct solution 


cotg 2 0 O * - 


1 _ 

«i 



1 ) • 


Note that here /3( = 1 — 
versa. 


-^ L ) is smaller than unity if 
M2 


<0 ( Pl >0), and vice 


With (D 3 - /3) <0 the reciprocal p - — leads to the monotonically con- 

w 

verging process 


p 




~ 1 + (/3 - D 3 ) Pj / 3 


where 



(32) 


is a suitable starting value. 

In all circumstances one (and only one) source level focus will be pro- 
duced by the third layer. The situation is similar to that depicted by the left 
diagram of Figure 3. 

The limiting process K 2 -* K t if applied to the expressions (7b) and (8) 
yield formulas that can be further rearranged by the use of the focal equation 
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(28) . The final results may be written as 


Xj = 2y, (K, + 1) J (2(^ + ^^^T^5) 

1 "fcT (1 + 6,q) 3 

W [(Kj + l) Yl ] 2 (i + q) 2 |_2/3 + ^ 1 + 6 iq : j|_2/3 + (1 + 6 lQ )^] ‘ 


B. Dg = 0 


This in fact is a two-layer case, as = Mi- Equation (6) is easily 

solved: 

2 

^ 1 + = D 3 ; q = cotg 2 0 O ' : = ^ 3 _ . 

The solution is real if the gradient in the third layer (p 3 ) is not negative and at 
the same time larger than the gradient in the layer below it (p 2 ) * There is no 
focus with 0 £ /u 3 < n 2 . The same rules govern the whole subdomain A < 0) , 
of which D 2 = 0 appears as a fitting terminal. However, in B (D 2 > 0) foci be- 
gin to appear in pairs if 0 s < p 2 , although there are still none if /u 3 is too 
small. This state of affairs persists through the entire B- and C -ranges and 
will be expected to be found with the value D 2 = 1 separating the two ranges. 

Expression (7b) and (8) assume the forms 


X F = 2 Y2 


I 


1 + K 9 

1 -k 2 


(Ps - 


i) 


OJ 



^"cT 

(l + q) 2 


If and p 2 are both zero. 
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c. Dj = i 


The gradient in the first layer is zero here, (M t = 0) . Three essentially 
different classes of velocity profiles (Fig. 7) obey the necessary condition (3) . 

The limiting process K t — 1 (Mi ^ 0, -* 0) transforms equation (6) 



FIGURE 7. TYPICAL VELOCITY PROFILES WHEN Mi = 0 


into 


v 3 + (y - 1) v + 7 (“T 2 - - 1) = 0 , 

M 3 


where 


y 


(K 2 + 1) & — . 

Yl 


(33) 


(34) 


For the iteration of equation (33) one uses better the unknown 


<P = v 3 = ( 1 + &>q) 3 / 2 


(35) 
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writing 


cp , - (1 - y) cp + y (i - ^ z ) 
^ v+l l> /x 3 


(36) 


The two classes of profiles at the left of Figure 7 are specializations of 
the classes B and C in Figure 3. As these do, they give rise to one single source 
level focus out of the third layer. However, as with pt 2 = 0, two iteration 
schemes must be set up to accommodate a coefficient in the equation, formerly 
(D 3 - /?), now (1 - y) . 

If the positive quantity y remains less than unity, the process (36) 
converges monotonically. A suitable initial value was found as 


<P i = 1 " 7 


Mi 

M 3 


(37) 


which is always positive under the circumstances envisaged. 
The rare event y = 1 gives the direct solution 
(i _J&)0 _ i 


q - cotg 2 d 0 " = 


Ms 


form 


With y > 1 the reciprocal of w , t = — , sets equation (36) into the 

<P 


V+1 r<‘-f> 

^3 


[ 1+ • 


(38) 


Monotonic convergence here is ensured; one may begin with 


ti = 


(1 - -^-) 3 
Ms 


(39) 
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As was expected, the right profile type in Figure 7 ( 0 ^ p 3 < p 2 ) does 
not always provide a third-layer focus. It is related to the B- and C-types of 
Figure 5. The critical quantity, there Q + R, here is 




(40) 


Focus formation is precluded if y §1, or, more generally, if S < 0. It cannot 
take place with too small a value of yu 3 . 


The special case S = 0 yields the explicit solution 
- 1 


q = 


1 - Y 
3 


^2 


With S > 0, when of necessity y < 1, two foci appear. One of them can 
be obtained after employing the approach (36) with the starting value (37) . The 
other one follows from equation (33) when it is given the form 


v + y , 
v Vita 


(ff - *) 


v+i 


1 - y 

The iteration may begin with 


(41) 


Vl 


JL 


[Mz 

1 - r 


“) 


(42) 


which is a positive quantity. 

The expressions (7b) and (8) for x and co must also be subjected to 

h 

the limiting process Kj -* 1. They then give 


2yi \Tq jjy + 2 + ^q]] 


1_ \ r q~ 

" = yi 2 ( 1 + q) 2 


A +. Aa. 


[y + 2 + ^q ] [y + 2 + 3 (^q] 
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Use has been made here of equation (33) for simplification. 


VH. A FOURTH SPECIAL TYPE: V 2 = V 0 


The three classes of profiles occurring here mark the boundary of the 
classes C and D on Figures 3 and 5 (where the case V 2 = V 0 , or D 2 = ^ had 
been excluded) . They are sketched in Figure 8 . 



FIGURE 8 . TYPICAL VELOCITY PROFILES WITH K 2 = 1 

Since K 2 = ^ =1, the inclination, 0 2 , of a ray when entering the third 
V ° 

layer * is equal to the inclination it had when leaving the sound source. With 
p 1 > 0 , a group of small elevation rays is sent back to earth from the first 
layer and cannot enter the third. All those that do enter it will also return to 
source level, since the gradient p 3 in the (infinitely extended) layer is not 
negative. 


With 62 = 1 - k 2 2 = 0 the u-equation ( 13) , derived from the original 
focal equation ( 6 ) , assumes the simple form 


0 2 is given by sin 0 2 = n/ 1 - K 2 2 cos 2 
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U(1 + D 2 - D 3 ) = D 2 • 


(43) 


If the parenthesis here is not negative, a positive solution for u exists. 
The condition may be written as 


> 0 . 

M; 3 Yl 

Since -^- L = 1 - -^ 2 -, it is equivalent to the requirement 
Yi 


(44a) 


m. < . y?. — 
~ y 2 - yi 


(44b) 


With this relation upheld the solutions u of equation ( 43) can be shown to be 
physically relevant; i. e. , they are larger or smaller than unity depending on 
whether > 0 or Sj < 0 * . 

Since condition (44b) is fulfilled by the profiles belonging to the two 
left types on Figure 8 one single ground level focus will always develop with these 
as it will with any profile from the subdomains C and D in Figure 3. 

The correspondence is not equally simple and complete when the third 
class of profiles is compared to the classes C and D on Figure 5. If m 3 is too 
small (so that condition (44b) is violated) , there will be no focus. An analogous 
behavior was found in the subdomain C; but in D, smallness of M3 does not pre- 
vent focus formation. On the other hand, if m 3 is sufficiently large, one single 

focus will appear with D 2 = , as it always will in subdomain D, while in C 

yi 

two foci will be generated by the third layer. 

The expressions to be derived for 9 0 " , x , cu from equation (43) are 

r 

relatively simple, and general results can be formulated as to the effects on 
them of m 3 - variations. In the following the height y 2 and the entire velocity 
profile below it will be considered as parametrically fixed. 


With 5 


1 


0, Vl = V 0 = v 2 , 


and q = 


Vfl_ 

M 3 y 2 


; see footnote on p. 26. 
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The quantity 


A = t _ Jf2. 22 JL 31 = i + J*L ll 
M3 Y2 M 3 Y2 


( 45 ) 


is not negatr, w-, on the strength of conditions (44b) or (44a) . Equation (42) 
can be solved for q to give 


q = cotg 2 0 O * = -j- - 1) . 


(46) 


One verifies without much labor that, in all three profile classes, a stronger 
third-layer gradient causes a steeper focal ray ( 0 O ' ,: increases) . 

Expression ( 7b) and ( 8) for the landing distance of the focal ray and 
for the intensity factor near the focus simplify to give 


2y 2 (K, +1) 


■R 


+ 6jq 


(47) 


w [y 2 (K t + l)] 2 


f L± g|Q 

, 1 + q 


\Tq" 


(48) 


For judging the dependence of these quantities on their argument q we 
consider the logarithmic differential quotients 


dx 
F 

dq 


2q(l + 6jq) 


> 0 


(49) 


1 = 2 6i - 1 J_ 

co dq (1 + djq) (1 + q) 2q 


(50) 


This distance x decreases with q, i.e. , with increasing angle of de- 

r 

parture ( 0 O ‘' ) . The steeper the focal ray, the nearer to the sound source it 
returns to ground level. 
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An equally unqualified statement cannot be obtained from the differential 
quotient (50) which can be positive, negative, or even zero. For the discussion 
we introduce o> from expression (48) and write 


do) 

dq 


1 

[y 2 (Kt + 1 )] 2 


i + frq 
2 'sTq ( 1 + q) 3 


[<5 iq (q+5)+ (i - 3q)] . 


(50a) 


The sign ambiguities in the last bracket are somewhat reduced by the 
restriction that cos 6 0 * must not be smaller than about 0. 95 ( see Introduction) , 
or that 


q = q « 9 • ( 51 ) 

The parenthesis ( 1 - 3q) is certainly negative then. 

A few words may be intercalated here on the significance of the restric- 
tion (51) regarding other quantities in the case V 2 = V 0 . First, although the 
landing distance decreases with q, it cannot decrease indefinitely without losing 
plausibility. To be sure, the focus would move into the sound source when 
q = 0, or 0 O * = 90° , but the supporting equations are then meaningless. Secondly, 
upper bounds are set up for the velocity gradients. From u 2 S 0 it follows that 



Um Vi 

Furthermore, since - = — — , 

Mi y 2 - yi 


M2 so. 0541 — ^ — 

y 2 - yi 


(46) and (45) , 
-■]" ■ 


Finally, on combining expressions 


M 3 = Mi 


h. 

y 2 


f — 4 

L 4i + <5jq 
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Evidently, the restrictions on Mi and p, 2 are relevant with profiles of the first 
type only (Mi > 0, M2 < 0) . The values of the three gradients encountered in 
practice rarely go beyond the above bounds, so that the restriction on 9 0 * sel- 
dom interferes with the interpretation of actual profiles. 


Returning now to the discussion of co let us first consider the left-hand 


class on Figure 8. Since Vi > V 0 here, the quantity 6 




and there- 


fore the differential quotient (50a) is negative. The intensity parameter co 
increases with 9 0 * : The larger M3, the steeper the focal ray, the closer the 
focus to the source, and the more serious is the focal situation. This is what 
might be called the regular sequence, also observed with the corresponding 
two-layer class arising when the middle layer is extended down to ground level. 


As an example, consider the profile where V 0 = V2 = 340 m/sec, Vi = 341 
m/sec, yi = 100 m, y 2 = 200 m. The upper bounds for Mi and - M2 are both 

0.184 — — , while, in fact, Mi = - A*? = 0. 01 — , well below the ceiling. Since 
sec ^ 1 sec ^ 

6i « - — — , the upper bound on M3 in « 0.188 . With 

x 1 u sec 


M3 

=0.02 : A 

= 1.25, 

m 3 

=0.04 : A 

= 1.125, 


61.2, 0* 
35.7, 0 O * 


7° 17', x p 
9°30', x p 


7. 84 km 


5. 39 km. 


When the focus moves from the first to the second position the average intensity 

near it increases by the factor ^ = 3. 33, or by 5. 2 db. 

wi 


Profiles with Mi > 0 confine low elevation rays to the ground layer. 

The first ray that can possibly escape it is characterized by \l 1 - K^cos^ 0 O = 9, 
orbycotg 2 0 o = - — . The gradient in the second layer being negative, that ray will 
traverse it and enter the third layer under the angle 6 = 0 O . It will be a focal ray 


(0 0 = dn' 1 ') if q = -~ satisfies the focal equation (46), i.e., if A = 00 and, 

therefore, M3 = 0. Because it has the lowest angle of departure of all the focal 
rays, it will be returned farther out and with lower intensity than any other ray. 
Indeed, by expression (47) and (48) , x = °° and co = 0. This clearly is a 


mathematical result devoid of practical significance. 
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With the two right-hand classes in Figure 8 (where Mi < 0, Sj > 0) 5 the 
ground layer cannot restrain any ray from going up and penetrating into the 
second layer. On the other hand, even a mathematical focal situation cannot 
develop if M3 = 0. By condition (44b) , the smallest value of M3 that still permits 
focus formation out of the third layer is given by 


Ml = Y z > i 

M 3 Y2 - Yl 

It defines a profile of the third class for which the zero-ray is focal ray, since 
A =0, q = 00 , 0 O = 0. This ray culminates at the boundary of the second and 
third layers and, by expressions ( 47) and ( 48) , is returned at the ( largest 
possible) distance 


X F = 2yz 


1 + Ki 
1 


with the intensity parameter 


O) 


°° . 


A very strong focus can be expected here. In passing, it will be observed that, 

d^ 

g 

if x and 2 „ _ „ * remain finite, the general expression (7) shows that 

F d»o “0 _ 

the intensity factor will always exceed all bounds if 0 O ‘ can and does approach 
zero. High noise levels at rather large distances may sometimes be explainable 
that way. 


On increasing M3 the focal distance will decrease; the intensity factor, 
contrary to its usual behavior, cannot but decrease, too, and does so at a fast 
pace, since the cotangent -function drops rapidly when its argument increases 
from zero on. This trend continues down to a minimum of a; which, by ex- 
pression (50a) , can be determined after solving the quadratic 


6 1 q 2 +(5<5 1 -3)q+l=0 
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Since the positive quantity 6 t is very small the solutions may be given as 


q « 


3 - 55 t r 

26 j [_ 





2 61 


(3 - 5^) 



or, more explicitly, as 


3 - 56, 

q « — - — L and q « 
°i 


1 

3-56! ' 


Of these, the second solution is too small and must be disregarded. The values 
of oj and A related to the first are found as 


1 

“ [y 2 (K t + 1)] 




9 


l/l 


] 


A 


m 


2 

4 - 56j 




2 * 


Rearranging relation ( 45) to give 


m = _J._ y?. - Yt , 

lh 1 “A y 2 

we see that enlarging M 3 results in larger values of A , and therefore of 6 0 . 

When M 3 has grown to equal m 2 > the dividing line between the two right-hand classes 

in Figure 8 has been reached, characterized by A = . If this value is still 

y 2 

smaller than A^.i.e. , if y 2 > 2yj, the intensity factor u> is still decreasing; if 

y 2 < 2y t , it will have passed its minimum and increase. Conceivably, the second 

layer's thickness, (y 2 - yi) , could be so large that no realistic value of M 3 will 
suffice to carry c 0 beyond the minimum. However, in ordinary circumstances 
values 

Y2 


36 



can be admitted. With such profiles the "regular" sequence obtains: c o becomes 
stronger with n 3 . 

Take, for instance, 

V 0 = V 2 = 340 m/sec, Vi = 339.32 m/sec, y t = 50 m, y 2 = 200 m . 

Figure 9 is based on these data. The focus abscissa, x , has been taken as the 

F 

independent variable instead of M 3 - The ration and the intensity level are 

M 2 

plotted versus x . 

F 

The second layer velocity gradient is the constant 

M 2 = ' | Mi , Mi = - 0. 0136 . 

The average intensity level near a focus, 


IL = 10 1og 10 , 

CO 

contains the arbitrary reference value . In Figure 9, co has been left un- 
determined, so that no absolute values are given on the db-scale. 

The foci at the left of x « 12. 2 km (M 3 = M 2 ) are created by profiles 
belonging to the middle class in Figure 8 . It is seen that, if-&. * 

1 . 5, the focus 
M 2 

will reside at x « 11 km, and the average intensity near it will be least alarm- 

r 

ing. It will, however, increase irrespective of whether M 3 is increased or de- 
creased. The "regular" (two-layer) sequence prevails at the left of the mini- 
mum. The curves are terminated at x =5 km, corresponding to M3 = 0. 0425 
( a value that seldom will be exceeded*; . The focal ray existing with this 


The sound propagation velocity would then increase by 4. 25 m/sec for every 
100 m up in the third layer. 
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FIGURE 9. GRADIENT RATIO REQUIRED TO LOCATE THE FOCUS AT x 
AVERAGE INTENSITY LEVEL THERE. K, = 1. 








































particular profile has the initial elevation 0 Q '' = 8° 20'. The maximum elevation, 

, is given by sin = n/ 1 - Kj 2 cos 2 0 O * ; since K* = ^ is but little differ- 

v ° , 

ent from unity, the angle 9{‘' is insignificantly larger than 0 O * and well within 
the range of validity for the 6. 

The foci at the right of x « 12. 2 km exist with third-class profiles 

(Fig. 8) . If the middle layer were extended to ground level, the corresponding 
two-layer atmosphere thus obtained could not support a focus at all. On 
decreasing ju 3 , the average intensity level grows very rapidly, and small dif- 
ferences in i u g (and x ) will affect it out of proportion. The focal rays here are 

all of very small initial elevation, the zero-ray ( 0 O : ' 5 = 0) arriving at the (largest 

possible) distance x = 12.64 km, corresponding to ^ = ^ 2 ' " = 0. 75 and 

F p 2 Y2 

co = °° . If ju 3 is decreased still further, a focal situation in the usual sense 
cannot develop. However, ray tubes of very small initial elevation will, what- 
ever the positive gradient p 3 might be, always arrive near x =12.64 km and 
enclose a very small frontal area causing a focus-like high intensity level. This 
is a general result: With profiles capable of returning the zero ray at a finite 
distance high levels will always be found there, although, in general, 9 0 = 0 
will not solve the focal equation. Machine computations in such cases show 
declining intensity with increasing small 0 o -values; however, since the first value 
is usually 0 O =0.2°, they do not exhibit the very large intensities to be expected 
if is even smaller. 


VIII. CONCLUSION 


The three-layer "planar" atmospheres can be classified with the use of 
the parameter 

D 2 = 1 - 

M2 

where the constant quantity denotes the sound speed variation per unit length 
climb in the k^ 1 layer. ^ 

Acoustical foci on the ground may develop if /jl 3 is positive. Adapted to 
the various atmospheric classes, computational schemes can be set up for finding 
the focus location and the average noise level in its immediate neighborhood. 
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If ju 3 > ju 2 one single focus will appear (provided the third layer is ex- 
tended enough to turn around the focal ray) . An analogous statement holds in a 
two-layered atmosphere with n 2 > 0 and ju 2 > Mi- Even if /u 3 < 0 , a second-layer 
focus can develop in those circumstances. 

In the realm 0 < M3 < M2 the existence criteria are less clear-cut. There 

can always be a third-layer focus if D 2 > ^ . For smaller values of D 2 , the 

yi 

value of M3 must attain a certain minimum before foci can appear; but then as 
many as two may develop with a given meteorology. These occurrences have no 
counterpart in the two-layered atmosphere where the equivalent condition 
0 < M2 < Mi does not permit any focus formation. 

A second-layer focus will move toward the sound source and the average 
intensity level near it will rise if at a given value of Mi the value of m 2 is increased. 
A third-layer focus does not always respond in the same manner to a variation 
of M3 at fixed values of Mi and M2- K may draw nearer to the sound source with 
diminishing intensity level. 

Systematic surveys of the behavior of third-layer foci are planned. These 
will be carried out with the relatively easy methods described in the present 
report. 
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